An Infinite Product for e'^ via 
Hypergeometric Formulas for Euler's Constant y 

Jonathan Sondow 

INTRODUCTION. Euler's constant, 7 = 0.5772156649. . . , is defined as the Umit 

7 = lim iH 1 1 Inn 

n^oo\ 2 n 

In this note we prove the infinite product formula for e'' = 1.7810724179. 
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where the nth factor is 
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Figure 1. The nth partial product a(n) of the infinite product for 

The path that led to (2) began with the proof [5] of the double integral formula for 
Euler's constant 
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an analog of double integrals [4] for ^(2) and ^(3), where ^(s) = I^„>i« ^ . A series 

in the proof (essentially the first series in (4) below), with term ratio a quotient of 
polynomials in n, gave the idea to write the series as a hypergeometric function (see [2, p. 
208]). A hypergeometric transformation then paved the way to (2). 

We give four proofs of (2). The first two are hypergeometric; the second uses the 
digamma function yr(t) = F' (t)/r(t) and is shorter. The last two proofs, which avoid 
hypergeometric s, depend on integral formulas for y : a classical one in the third proof, 
and (3) in the fourth. Using the classical integral, we end the note with a quick proof of 
(3). 

PROOF 1. We follow the latter part of the path described. Using ln((n + l)/n) ^ to 
replace Inn by ln(n + 1) in (1), we convert the limit to the series 
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Interchanging summation and integration (justified by uniform convergence, the 
integrand being bounded by n~^{n + 1)"^), we write the formula as 



7= ) 9 dt=\ > dt 
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in order to replace the series by a hypergeometric function [2, p. 205], [3, p. 101] 
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, s-t s(s + l)-t(t + l) 
= 1 + — + ^ - + ■ 

u • V u(u + 1) • v(v + 1) 
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obtaining 
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The transfonnation 
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(valid if 5 > and v — t>0, and proved by applying the identity F\' M = FV' '\, 

then Thomae's transformation [3, pp. 104-105, 111], then the identity) gives the simpler 
expression 
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Expanding by (5), and simplifying terms, we get the key formula 
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Integrating termwise (the nth term being bounded by (n + 1) ^) by the method of partial 
fractions, we arrive at the series 



°° 1 n 



\n(k + l). 



and exponentiation yields the desired product (2). 



PROOF 2. Using values of the y/ function and a series for y/' (see [1, PP- 15, 16, 22, 
equations (4), (8), (9), resp.]), we shorten the first proof. From (1) and the formulas 
= ~7 ¥(n + 1) = 1 + 1/ 2 H h 1/ n - 7 , we deduce that 



f [r(0-Jjrfr = r. 

Employing (5) and (6), we transform the series i/a' (t) = £„>o(? + «)~^ into the function 
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Using (5) (but not (6)), we verify the identity 
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and (7) follows. The rest of the proof is the same. • 

Remark. One can accelerate the convergence of such hypergeometric representations by 
computing some terms in (1). For example, 
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dt. 



PROOF 3. We reduce the key formula (8) to the classical formula for Euler's constant 

r = \li7^ + ^]du. (9) 



(One way to derive (9) is by inserting the value y/(l) = -y into Gauss's integral for the \\f 
function [6, p. 246, Example 1].) Denoting integral (8) by /, we substitute Euler's beta 
integral 
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(evaluated by n-fold integration by parts), then replace r by r + 1, obtaining 
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For fixed r > , a little calculus shows that 
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guaranteeing uniform convergence of the series 
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on the interval 0<m<1. This permits interchanging summation and integration with 
respect to u, which gives 



i=rCui^-i]dudt. 
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(10) 



Reversing the order of integration (the integrand being non-negative), we integrate with 
respect to t, and arrive at integral (9). This proves (8), and (2) follows as before. • 
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PROOF 4. We show that integrals (8) and (3) are equal. Repeating the previous proof 
through (10), we substitute the integral 



reverse the order of integration from dvdudt to dtdvdu (non-negativity again), and 
integrate with respect to t, obtaining 



Since the change of variables v = 1- x,u = xy in integral (3) transforms it into integral 
(11), we infer (8), and (2) follows as usual. • 

Remark. Integration with respect to v in integral (11) produces integral (9) for 7. 
Combined with the equality of integrals (3) and (11) by change of variables, this gives a 
short proof of (3). (The proof in [5] is longer, but uses only the definition of Euler's 
constant.) 
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